A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

11(11.07.1; 11.19.6); 12(12.12.1); 13(13.18.1) 



ASTRONOMY 

AND 
ASTROPHYSICS 

February 1, 2008 



The Phoenix radio survey: the angular correlation function 

A. Georgakakis^'^, B. Mobasher^, L. Cram"^, A. Hopkins'*, and M. Rowan-Robinson^ 

^ School of Physics and Astronomy, University of Birmingham, Edgbaston, B15 2TT, UK 

Astrophysics Group, Blackett Laboratory, Imperial College, Prince Consort Rd , London SW7 2BZ, UK 
^ Astrophysics Department, School of Physics, University of Sydney, NSW, Australia 2006 
^\ ' * Department of Physics and Astronomy, University of Pittsburgh, 3941 O'Hara Street, PA 15260, USA 

as: 

CT^ • the date of receipt and acceptance should be inserted later 



o 
O 
oo 



> 

oo 

m 

o 

ON 
ON 

I 

o 



% 



Abstract. Using the Phoenix radio survey, a homoge- 
neous survey selected at 1.4 GHz and covering an area 
of «3deg^, we analyse the clustering of the sub-mJy ra- 
dio population using angular correlation function analy- 
sis. Extensive simulations are carried out to investigate 
the significance of the estimated angular correlation am- 
plitudes. Our analysis show that for the 5i.4 > 0.5 mJy 
sub-samples the radio source distribution is anisotropic 
at the 2(7 significance level. Additionally, we estimate up- 
per limits for the angular correlation amplitudes that, de- 
spite the large uncertainties, are in good agreement with 
the amplitude estimates for sources brighter than 1 mJy, 
detected in the FIRST radio survey (Cress et al. 1997). 
Adopting a radio luminosity function and assuming an 
evolving spatial correlation function of the form ^(r) = 
('''/ro)~''(l -I- 2)~^'^+^-', with the evolution parametrised by 
e, we find an upper limit for the angular correlation length 
To ~ 9h~^ Mpc for 5*1.4 > 0.5 mJy and 7=2.1. This agrees 
well with the value Tq ~ 6 — Sh^^Mpc estimated from 
the FIRST radio survey for sources brighter than 1 mJy. 
Additionally, we quantify the characteristics, in terms of 
areal coverage and limiting flux density, of future deep 
radio surveys to yield a significant correlation amplitude 
detection and to explore possible changes of the correla- 
tion amplitude with flux density. 

Key words: Galaxies: general - Galaxies: structure - 
Cosmology: large scale structure of Universe - Radio con- 
tinuum: galaxies 



1. Introduction 

The spatial and angular correlation functions, ^(r) and 
'w{9) respectively, have been extensively used to quantify 
the clustering of galaxies. The former has the important 
property of being the Fourier conjugate of the power spec- 
trum, P{k), of the galaxy distribution (e.g. Peacock & 
Nicholson 1991). The latter is the two dimensional projec- 
tion of ^(r) on the sky and is related to it by the Limber 
equation (Limber 1953; Phillipps et al. 1978). Although 
w{9) conveys less information on galaxy clustering than 



^(r) due to its integral nature, nevertheless it is a powerful 
probe of the large scale distribution of galaxies, especially 
at faint magnitudes where redshift information is difficult 
to obtain. 

Wide-area surveys at optical wavelengths have pro- 
vided information on the large-scale structure at redshifts 
z « 0.1 (Loveday et al. 1995; Maddox et al. 1990b). 
Additionally, deep optical and near-infrared samples over 
smaller areas (Villumsen et al. 1997; Hudon et al. 1996; 
Roche et al. 1998a, b; Carlberg et al. 1997) give use- 
ful information on the large-scale distribution of galax- 
ies at higher redshifts (z <1.0). There is a general con- 
sensus that at optical wavelengths, w{9) is a power law 
of the form w{e) = A,,e~('^-'^\ with 7 w 1.8 (Mad- 
dox et al. 1990b). The amplitude, A^, decreases with in- 
creasing depth of the survey, while there is also evidence 
that the slope, 7, flattens at faint magnitudes (Infante & 
Pritchet 1995). The spatial correlation function has the 
form ^(r) = {r/ro)~^, where the correlation length, Tq, 
is found to be r^ = 5.4h-iMpc (Davis & Peebles 1983). 
However, it has been demonstrated that early-type galax- 
ies are more clustered (higher values of Tq) than late types 
and that lower luminosity galaxies are a factor of ~ 2 
less clustered than their brighter counterparts of the same 
Hubble type (Loveday et al. 1995). There is also evidence 
that To decreases with increasing depth of the survey (In- 
fante & Pritchet 1995), due to the presence of a population 
of weakly clustered faint galaxies. 

Radio surveys, unlike the optical ones, are not affected 
by galactic dust extinction and mainly consist of high 
redshift objects (z « 1). Therefore, they sample much 
larger volumes than optically selected samples, albeit with 
sparser coverage, providing the opportunity to study the 
clustering of matter at much larger physical scales. How- 
ever, the broad luminosity function of extragalactic radio 
sources (Condon 1984; Dunlop & Peacock 1990) implies 
that a flux density-limited sample spans a wide range of 
redshifts. The projection on the sky of all detected objects 
results in a distribution close to random, smearing out 
any information on the large-scale structure. This is sup- 
ported by the fact that no signal has been detected in the 
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angular (2-D) correlation analysis of bright radio galaxies 
(Webster 1976; Masson 1979). However, Shaver & Pierre 
(1989) found evidence of anisotropic distribution of ra- 
dio sources towards the supergalactic plane extending out 
to at least z w0.02. More recently, a clustering length of 
To = 11 h^^ Mpc was estimated at 1.4 GHz from the spa- 
tial (3-D) correlation analysis of radio sources having flux 
densities > 0.5 Jy and redshifts in the range 0.01 < z < 0.1 
(Peacock & Nicholson 1991). 

Cress et al. (1996) have suggested that at lower flux 
densities, projection effects may become less significant 
and hence, a 2-D correlation analysis can be applied 
successfully to investigate the clustering of radio galax- 
ies. Indeed, a non-zero amplitude has been found for 
the angular correlation function of radio sources with 
Si.sbGHz > 50 mjy in both the Green Bank (northern 
hemisphere) and the Parkes-MIT-NRAO (southern hemi- 
sphere) 4.85 GHz sky surveys (Loan, Wall & Lahav 1997; 
Kooiman, Burns & Klypin 1995). More recently, using the 
FIRST radio survey (1.4 GHz; Becker, White & Helfand 
1995), with a uniform sensitivity of ImJy over an area 
of 1500 deg^, a non-zero and clearly significant amplitude 
for the angular correlation function is estimated (Cress et 
al. 1996). Adopting the radio luminosity functions (RLF) 
determined independently by Condon (1984) and Dunlop 
& Peacock (1990), Cress et al. (1997) inferred a clustering 
length of ro=6-8h^^ Mpc for Sia >lmJy radio sources. 

The recent large-area radio surveys described above, 
have provided information on the two-dimensional pro- 
jected distribution of relatively bright (51.4 >1 mJy) radio 
sources, dominated by bright ellipticals and AGNs. How- 
ever, there is still limited information on the clustering 
properties of the faint (sub-mJy) radio population. At flux 
densities below few mJy there is evidence for the appear- 
ance of a new population of faint radio sources, likely to 
comprise a large fraction of starbursts (Benn et al. 1993; 
Georgakakis et al. 1999). In particular, the radio lumi- 
nosity function models developed by Condon (1984) pre- 
dict that the surface density of starbursts increases from 
10% of the radio population at S'1.4 « 1 mJy to over 30% 
at 51.4 ~ 0.4 mJy. Moreover, studies of the spatial dis- 
tribution of starbursts, selected at optical and infrared 
wavelengths, show that these objects are expected to have 
different clustering properties (Davis & Geller 1976; Gio- 
vanelli et al. 1986; Saunders, Rowan-Robinson & Lawrence 
1992; Loveday et al. 1995) compared to those of early-type 
galaxies, with which bright radio sources are often associ- 
ated. Therefore, study of the angular correlation function 
at sub-mJy flux density levels has the potential to reveal 
differences in the clustering properties of the faint radio 
population. 

Recently, Benn & Wall (1995) demonstrated that the 
isotropy (or anisotropy) in the radio source counts in dif- 
ferent fields of similar geometry and sensitivity can be used 
to set limits on the scale of the largest cellular structures in 
the Universe. However, they argued that at sub-mJy and 



/ijy flux densities, such a study is hampered by the small 
number of comparable surveys and by the small number 
of sources detected in each field. This is also supported 
by Windhorst et al. (1990), who reviewed the field-to- 
field variations in the radio source counts of small area 
radio surveys (< Ideg^). They concluded that the differ- 
ences, although exceeding the random distribution expec- 
tation, are due to statistical fluctuations, rather than of 
cosmic nature. Nevertheless, a non-uniform angular cor- 
relation function was estimated by Oort (1987a) for the 
radio sources detected in the deep {Sia > 0.1 mJy) Lynx 
fields (Oort 1987b). More recently, Richards et al. (1999) 
also reported the detection of clustering signal for radio 
sources brighter than 40 /iJy, detected within a 40 arcmin 
diameter radio survey (1.4 GHz) centred on the HDF. 

The Phoenix radio (1.4 GHz) survey, covering an area 
of 3.14 deg^, larger than any other survey at a similar flux 
density limit (S'1.4 = 0.4 mJy) and reaching surface densi- 
ties « 6.6x 10^ sources sr^^, provides a unique opportunity 
to study the clustering properties of the faint radio popu- 
lation. In this paper the angular correlation function of the 
radio sources detected in the Phoenix field is estimated. 
Section 2 gives a brief description of the observations. The 
method for calculating w{9) is outlined in section 3. In 
section 4 we estimate the correlation function of the ra- 
dio sample. Section 5 presents the simulations carried out 
to investigate the significance of our results, while in sec- 
tion 6 the correlation function amplitudes are compared 
with 3-D clustering models. The results from the radio 
correlation analysis are discussed in section 7. Finally, we 
summarise our conclusions in section 8. 



2. Radio Observations 

The Phoenix survey has been performed at 1.4 GHz us- 
ing the 6A configuration of the Australia Telescope Com- 
pact Array (ATCA). It is a mosaic of 30 pointing cen- 
tres covering a 2° diameter field centred at RA(2000) = 
O^h i^m ^2n6 ; Dec.(2000) ^ -45° 44' 8'.'0. The final 
synthesised beam size is «8arcsec. 

Details of the observation, data reduction, source de- 
tection and catalogue generation are presented in Hopkins 
et al. (1998). A source is included in the sample if its peak 
flux density is Aa above the local RMS and if it survives vi- 
sual inspection. There are two sources of incompleteness 
in the catalogue, as with any peak flux density limited 
sample. The first is a result of the attenuation of the pri- 
mary beam sensitivity away from a field centre. This effect 
has been minimised in the Phoenix survey by the mosaic- 
ing strategy used. The second is the fact that extended 
objects with a total flux density above the survey limit 
can be missed by a source detection algorithm which ini- 
tially detects candidates based on their peak flux densities. 
These effects have been described in Hopkins et al. (1998). 
The catalogue is estimated to be 50% complete to 0.3 mJy, 
and 80% complete to 0.4 mJy Therefore to minimise the 
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effect of incompleteness when performing the correlation 
analysis, we restrict ourselves to a subsample containing 
all sources brighter than 0.4 mJy. 

3. Two point correlation function 

The two-point angular correlation function, w{9), is de- 
fined as the joint probability 5P of finding sources within 
the solid angle elements Sfli and 5512, separated by an 
angle 9, in the form 



sp = N^ {1 + w{9))snidn2, 



(1) 



where N is the mean surface density of galaxies. For a 
random distribution of sources w{9)—0. Therefore, the an- 
gular correlation function provides a measure of galaxy 
density excess over that expected for a random distribu- 
tion. Various methods for estimating w{6) have been in- 
troduced, as discussed by Infante (1994). In the present 
study, a source is taken as the 'centre' and the number 
of pairs within annular rings is counted. To account for 
the edge effects, Monte Carlo techniques are used by plac- 
ing random points within the area of the survey. The most 
commonly used estimators in this procedure have the form 



DD 

w = „„ 
RR 

DD 

w = „ „ — 1, 
DR 



1, or 



(2) 
(3) 



where DD and DR are respectively the number of sources 
and random points computed at separations 9 and 9 + d9 
from a given galaxy. Similarly, RR is the number of ran- 
dom points within the angular interval 6 to + d6 from 
a given random point. Infante (1994) and Hcwctt (1982) 
emphasised the importance of correcting for any spurious 
cross correlation between the random and the galaxy cat- 
alogues when using the above estimators, arising partly 
because of the anisotropic and inhomogeneous galaxy dis- 
tribution relative to the field boundaries. The corrected 
form of the estimator is then 



DD 

Ur 



RD 

'rr' 



(4) 



where RD is the number of random-data pairs, taking the 
random points as centres (Infante 1994). Landy & Szalay 
(1993, hereafter LS) introduced the estimator 



DD - 2DR + RR 
RR ■ 



(5) 



The advantages of this estimator is that it minimises the 
variance of w{9) and also reduces the edge effects from 
which both the estimators (^) and ^ suffer. Moreover, 
it has similar properties (Landy & Szalay 1993) to that 
introduced by Hamilton (1993) 



DDRR 
'"=DRDR-^- 



(6) 



In the present study the correlation function has been cal- 
culated using estimators (|^), (||) and (^). The results agree 
within the respective uncertainties. Therefore, only the re- 
sults from the LS estimator are presented. 

For a given radio galaxy sample, a total of 100 random 
catalogues were generated, each having the same number 
of points as the original data set. The random sets were 
cross-correlated with the galaxy catalogue, giving an av- 
erage value for w{9) at each angular separation. In pro- 
ducing random sets of points, we take into account vari- 
ations in sensitivity, which might affect the correlation 
function estimate. The flux density threshold for detec- 
tion depends on the local RMS noise which varies across 
the survey area. Since the random fields are expected to 
have the same sources of bias as the data (i.e. the simu- 
lated catalogues must have the same selection effects as 
the real catalogue), the RMS noise map, with a resolution 
of ~10arcsec, is used to discard simulated points in noisy 
areas. This is accomplished, to the first approximation, by 
assigning a flux density to random points using the Wind- 
horst et al. (1985) logTV-logS* distribution, following the 
method of Cress et al. (1996). If the flux density of the 
random point is less than 4 times the local RMS noise, 
the point is excluded from the random data set. However, 
producing random points with a uniform distribution over 
the observed area does not change our final result. 

The uncertainty in w{9) is determined using both Pois- 
son statistics and 50 bootstrap resamples of the data 
(Ling, Frenk and Barrow 1986). For the latter method, 
simulated data-sets were generated by sampling N points 
with replacement from the true data-set of N points. The 
correlation function is then calculated for each of the boot- 
strap samples, following the same procedure as that with 
the original data-set. The standard deviation around the 
mean is then used to estimate the uncertainty in the cor- 
relation function. 

Although the LS estimator used here is shown to have 
Poissonian variance for uncorrelated points (Landy & Sza- 
lay 1993), it docs not necessarily behave this way for cor- 
related data. The bootstrap method is beheved to give 
a more representative estimate of the uncertainty associ- 
ated with w{9). However, Fisher et al. (1994) carried out 
a detailed study of the biases affecting the bootstrap er- 
rors (e.g. cosmic variance, sparse sampling by galaxies of 
the underlying density distribution) and concluded that 
overall the bootstrap uncertainties overestimate the true 
errors. Nevertheless, the bootstrap resampling is a general 
method for assessing the accuracy of the angular correla- 
tion function estimator and it will be used here to calculate 
the formal errors. The uncertainty associated with w{9), 
estimated by the bootstrap resampling technique, is found 
to be about three times larger than the Poisson estimate. 

There is expected be a significant number of physical 
double sources in a radio survey with the angular res- 
olution of the Phoenix survey (« 8"). In any study of 
the clustering of radio galaxies via correlation analysis. 
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Flux density (mjy) Number of Sources 

>0.4 634 

>0.5 529 

>0.6 454 

>0.7 391 

>0.9 316 

[0.4,0.9] 318 

Table 1. Number of sources in each flux density-limited 
subsample 



these should not be considered as two sources since both 
components are formed in the same galaxy. To identify 
groups of sources that are likely to be sub-components 
of a single source we have employed a percolation tech- 
nique where all sources within a given radius are replaced 
by a single source at an appropriate 'centroid' (Cress et 
al. 1996; Maghocchetti et al. 1998). Following the method 
developed by Maghocchetti et al. (1998), we vary the link- 
length in the percolation procedure according to the flux of 
each source. In that way, bright sources are combined, even 
if their angular separation is large, whereas faint sources 
are left as single objects. This technique is based on the 
9 — S relation for radio sources, which has been shown to 
follow e(x^/S (Oort 1987a). 

To define the relation between flux density and link- 
length, the angular separation of double sources versus 
their total flux is plotted in Figure |l|, out to a separa- 
tion of 180arcsec. Visual inspection has conflrmed that 
the pairs on the left of Figure H are predominantly sub- 
components of a single source. This is based on an as- 
sessment of the appearance of the object, including the 
disposition of the sources, the nature of any bridging ra- 
diation, and the appearance of source edges. Accordingly, 
we set the maximum link-length to be 



^link 



20 X ./Ft 



total 5 



(7) 



where Ftotai is the total flux of each group. This is shown 
by the dashed line in Figure |] and effectively removes the 
majority of visually identified doubles. Moreover, one can 
apply an additional criterion to identify genuine doubles, 
based on the relative flux densities of the sub-components 
(Maghocchetti et al. 1998). This is because lobes of a single 
radio source are expected to have correlated flux densities. 
Here, the groups of sources identified by the percolation 
technique are combined only if their fluxes differ by a fac- 
tor of less than 4. This procedure was repeated until no 
new groups were found. The final catalogue consists of 908 
objects to the limit of 0.1 mJy, with a total of 30 groups 
of sources being identified and replaced. 

The source counts of the sample, normalised to the 
Euclidean slope, are plotted in Figure along with the ra- 
dio counts at 1.4 GHz taken from Windhorst et al. (1993). 
There is a drop in the number counts below 0.5 mJy, as our 
sample is affected by incompleteness. To quantify this, we 
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Fig. 1. Angular separation against total flux density of 
double sources. The dashed line represents the maximum 
link-length, for a given flux density, used in the percolation 
technique. 

firstly fit a straight line to the source counts of Windhorst 
et al. (1993) for flux densities fainter that 5mJy (continu- 
ous line in Figure 0). We then compare our number counts 
at a given flux density bin with those predicted by the fit- 
ted line. We conclude that our sample is « 80% complete 
to 0.4 mJy. This is in agreement with the correction fac- 
tors for incompleteness independently derived by Hopkins 
et al. (1998) accounting for both resolution effects and 
the attenuation of the beam away from the field centre. 
Therefore, to minimise the effect of incompleteness of the 
radio catalogue, when performing the correlation analysis, 
we restrict ourselves to a subsample containing all sources 
brighter than 0.4 mJy. 

Finally, before fitting a power law to w{9), we take 
into account a bias arising from the finite boundary of 
the sample. Since the angular correlation function is cal- 
culated within a region of solid angle il, the background 
projected density of sources, at a given flux density limit 
So, is effectively NsiS > So)/fl (where N ^ (S > So) is the 
number of detected objects brighter than 5*0). However, 
this is an overestimation of the true underlying mean sur- 
face density, because of the positive correlation between 
galaxies in small separations, balanced by negative values 
of wlO) at larger separations. This bias, known as the in- 
tegral constraint, has the effect of reducing the amplitude 
of the correlation function by 



uJn 



^ J J w{e)dnidn2, 



(8) 



where il is the solid angle of the survey area. The angu- 
lar constraint is estimated using Monte Carlo integration. 
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Fig. 2. Circles: source counts at 1.4 GHz for the Phoenix 
Survey; Stars: source counts from Windhorst et al. (1993); 
Line: best fit to the Windhorst et al. (1993) source counts 
for flux densities less than 5 mJy. 

assuming the correlation function to follow a power law 
w{e) =A^x 0-^ with (5=0.8 or 1.1. We find uJn=lA6A^ 
and ojn—1.89Ai^ for (5=0.8 and (5=1.1 respectively, where 
9 is measured in degrees. 

4. Correlation Function of the Phoenix radio 
survey 

We have extracted from the original catalogue four flux 
limited samples with flux density cutoffs at 0.4, 0.5, 0.6 
and 0.7 mJy. Furthermore, to investigate changes in the 
clustering properties of radio sources with fiux density, 
two independent sub-samples were considered, having the 
same number of objects and fiux densities in the range 
0.4 < 51.4 < 0.9 mJy and Sia >0.9mJy Table lists 
the number of sources in each sub-sample. The correlation 
function w{9) is calculated for each sub-sample for angu- 
lar separations ranging from 0.03 deg to 1.3 deg, within 9 
equally separated logarithmic bins. The results are shown 
in Figures and 0. 

The amplitude of the correlation function is estimated 
by fitting the following function to the observations 



w{e)^Ay, xe-' -LUniS) 



(9) 



where the integral constraint is evaluated for the two dif- 
ferent values of i5 adopted in previous studies; (5=0.8 (Pea- 
cock & Nicholson 1991; Loan, Wall k Lahav 1997) and 
(5=1.1 (Cress et al. 1996). The calculated amplitudes and 
the Poissonian uncertainties are listed in Table y. Assum- 
ing Poisson errors, these amplitudes are non-zero at the 2(7 



Density (mJy 


A':='-'> 




>0.4 


O.OOliO.OOl 


0.003±0.003 


>0.5 


0.003±0.001 


0.006±0.003 


>0.6 


0.004±0.002 


0.008±0.004 


>0.7 


0.005±0.002 


0.010±0.004 


>0.9 


0.004±0.002 


0.009±0.005 


[0.4,0.9] 


0.003±0.002 


0.007±0.005 



Table 2. Angular correlation function amplitudes and the 
associated Poisson errors. 



confidence level for all the sub-samples, except from those 
with S'i.4 > 0.4 mJy and 0.4 < S'1.4 < 0.9 mJy. However, 
since the formal bootstrap errors are « 3 times larger than 
the Poisson expectation, the amplitudes derived here are 
treated as upper limits. 

Additionally, there are two effects that might further 
increase the correlation amplitude errors: (i) the inter- 
dependence of the w{6) measurements at different angles 
6 and (ii) the dependence of the uncertainty estimates on 
higher order correlations (skewness and kurtosis) that are 
ignored in this study. However, Mo et al. (1992) demon- 
strated that the underestimation of the correlation ampli- 
tude uncertainty, (5^^, , by these two effects is compensated 
by the overestimation of the ^(6*) errors by the bootstrap 
resampling method. 

Furthermore, because of the small solid angle of the 
Phoenix field, cosmic variance may affect the estimated 
correlation amplitudes. The effect of cosmic variance is 
studied in detail in the next section. 



5. Simulations 

Because of the relatively small Phoenix field size, one 
should be cautious about the interpretation of the cor- 
relation amplitude upper limits, since field-to-field fluctu- 
ations may affect our results. For example Oort (1987a) 
carried out the angular correlation function analysis using 
sub-mjy radio surveys smaller than Phoenix and found 
evidence for anisotropic source distribution in only some 
of them. 

Firstly, to test the significance of our results we con- 
struct catalogues of randomly distributed points and cal- 
culate the correlation function from a region with the 
same geometry as the Phoenix field. In all trials, the es- 
timated amplitudes are found to be consistent with zero 
within the Poisson standard deviation, crpoisson (Figure 
W. Therefore, the amplitudes derived here (non-zero at 
the 2 a Poisson significance level) , imply a non- uniform dis- 
tribution of the faint radio population, albeit at the 2a 
confidence level. The exceptions are the S'1.4 > 0.4 mJy 
and 0.4 < S'1.4 < 0.9 mJy sub-samples. This might indi- 
cate a low correlation amplitude at faint flux densities (e.g. 
Au, < 3 X 10"^, 7 = 1.8, Si. 4 > 0.4 mJy), that cannot be 
detected by the small size of the Phoenix field. However, 



A. Georgakakis et al.: Correlation function of the Phoenix radio survey 




-O.B -0.4 

log(e) (deg) 




-O.B -0.4 

log(e) (deg) 



_ Flux Density Cutoff : 0.4mJy 



_ Flux Density Cutoff : O.SniJy 



-J— ^- 



-0.8 -0.4 

log(e) (deg) 




-0.8 -0.4 

log(e) (deg) 



Fig. 3. Angular correlation function w{9) of the Phoenix radio survey for subsaniples with flux density cutoffs at 0.4, 
0.5, 0.6 and 0.7 mJy. The error bars shown are Poisson estimates. The solid and dashed lines represent the fits to the 
data assuming (5=0.8 and 6=1.1 respectively. Points in the range 3' < 9 < 20' are used in the fitting algorithm. 



one should be cautious about this interpretation, since at 
0.4 mJy the sample is also likely to be affected by incom- 
pleteness. 

To further investigate the sensitivity of our results 
to cosmic variance, artificial galaxy catalogues are con- 
structed with a built-in correlation function of the form 
Aiij 9^ . Then, we attempt to recover the correlation am- 
plitude from a field having the same geometry as the 
Phoenix survey. 

The algorithm used to generate simulated catalogues 
with projected hierarchical power-law clustering is de- 
scribed by Infante (1994) and is similar to that employed 
by Soneira & Peebles (1978). Firstly, two points separated 
by an angle 9i are randomly placed within a 50x50 de- 
grees area. Each of these points serves as the center of a 
new pair of randomly oriented points with angular sep- 
arations ^2 ~ ^i/-^; where A is a real number. The four 
points generated at the previous step are the new centers 
for pairs of points separated by 9^ = 92/X. Therefore the 



Lth step produces 2^ ^ pairs with angular distance 6l 
between the points of the pair. 

Adopting a value A ~ 9i/9i^i = 1.8 for the ratio of 
separations between successive levels, produces a power- 
law distribution of points with exponent 6 « 0.8 (Soneira 
& Peebles 1978). The angular distance between points at 
the first level is set 01 — 1.5 deg. Additionally, the number 
of hierarchical levels, L, is drawn from a uniform proba- 



bility distribution in the interval [L 



minimum and maximum number of levels are such that 
the smallest angular separation of the points of a pair is 
less than «1 arcmin. However, these parameters produce 
too many points at each level, resulting in large angular 
correlation amplitudes, Ay^, compared to those found in 
bright radio samples (Loan, Wall & Lahav 1997; Cress et 
al. 1996). Therefore, each point is assigned a fractional sur- 
vival probability, /. This step does not change the form of 
the built-in angular correlation function (Bernstein 1994), 
but allows tuning of A^- Values of / between 0.01 and 
0.05 produce amplitudes in the range 1 x 10"'^ to 20 x 10~^ 
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Fig. 4. Angular correlation function for sources with flux densities 5i.4 >0.9mJy and 0.4 < Si, 4 < 0.9niJy. The error 
bars shown are Poisson estimates. The solid and dashed lines represent the fit to the data assuming (5=0.8 and (5=1.1 
respectively. Points in the range 3' < 9 < 20' are used in the fitting algorithm. 







' 1 ' 


1 1 1 1 1 








Trial n 








A,-(10.0^0,6)xl0 3_ 





- 




- 










LO 





>. 


x* : 














- 


1 


1 ,1 



Trial D 
A^-(4,3±0,5)xl0-: 



-1.2 -1 -O.H -0.6 -0.4 

log(0) (deg) 



-1.2 



-1 -0.8 -0.6 -0.4 
log(e) (deg) 



Fig. 6. The angular correlation function for the simulations B and D calculated from an area ~ 10 times larger than 
that of the Phoenix field. The error bars are Poisson estimates. 



respectively. Here we only consider amplitudes in the in- 
terval 1 X 10^^ < Au, < 10 X 10^^, similar to those found 
in brighter radio samples (Cress et al. 1996). The param- 
eters used for four such simulations (trials A-D) are listed 
in Table I . The estimated 'w{9) for trials B, D are plotted 
in Figure ^. 

The procedure described above is repeated until the 
surface density of points within the 50x50 degrees area 
exceeds the surface density of sources of the present radio 
sample, at a given flux density cutoff. The angular corre- 
lation function is then calculated from the central 36 deg^ 



region, to ensure a uniform surface density of points. The 
amplitude of the correlation function is estimated by fit- 
ting the following function to the observations 



w{0) = A^"' X 



-0.8 



^n, 



(10) 



where the integral constraint, ujn, is 0.54 A^'™ for 6 = 
—0.8. The fitting is performed for angular separations 
between 3 to 50 arcmin. Additionally, for a given mock 
catalogue, the angular correlation function is also calcu- 
lated for nine fields, distributed within the central region 
of the 50x50 degree area, each having the geometry of 
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Fig. 5. The angular correlation function for a random 
sample of points. The error bars are Poisson estimates. 




Fig. 

0- 



7. Cumulative distribution of A defined in equation 



Trial A 

A 1.8 

B 1.8 

C 1.8 

D 1.8 



J-'min 

6 



,f 

0.03 
0.04 
0.01 
0.02 



density (deg ^) 
145 
170 
170 
200 



Table 3. Parameters used to construct simulated cata- 
logues. The densities of 145, 170 and 200 sources per deg^, 
correspond to the surface densities of the present sample 
at flux density cutoffs of 0.6, 0.5 and 0.4 mJy respectively. 



the Phoenix survey. The amplitude, A^***, and the Poisso- 
nian uncertainty, apoisson , are calculated (see section |3|) 
for each of these smaller fields and then compared to A^"^. 
The difference. A, between A'^* and A"^™ normalised to 
the Poisson standard deviation is defined 



A 



\At 



A':. 



<7Po 



(11) 



For a given survival probability, /, and surface density of 
objects (corresponding to a flux density cutoff) six artifi- 
cial catalogues are generated. The cumulative distribution 
of A values (Figured), shows that in w 70% of the trials, 
A^^* lies within « "^apoisson from Af„™. For a Gaussian 
distribution, this is the definition of one standard devia- 
tion. Therefore, the simulations carried out here indicate 
that the angular correlation amplitudes determined from 
a radio survey with similar characteristics (i.e. geometry, 
surface density of objects) to those of the Phoenix sur- 
vey, are representative of the faint radio population ampli- 
tudes within two Poisson standard deviations. This error is 
smaller than that estimated by the bootstrap resampling 
technique (~ Sapoisson'i see section H), suggesting that 
cosmic variance is not significantly affecting the present 
study. However, it is likely that the simulations presented 
here do not correctly model the higher order correlations 
produced by gravitational instabilities, on which the w{9) 
errors depend. Nevertheless, the Poisson errors estimated 
here are independent of higher order correlations and are 
only meaningful when compared relative to the Poisson 
uncertainties of the real data-set. In the rest of this study 
we retain the bootstrap uncertainties as the formal errors, 
since this method has been widely used to assess the ac- 
curacy of w{9). 

6. The Model 

Let us assume that the evolution of clustering of galaxies 
with redshift is parametrised by e such that their spatial 
correlation function is (Phillipps et al. 1978) 






(1 + 4 



<3+e) 



(12) 



where r is the proper length and Tq the correlation length 
for z = 0, corresponding to the scale where the clustering 
becomes non-linear. The clustering evolution parameter 
has the property that if e > the physical length of a 
typical cluster reduces with time, resulting in a growth of 
clustering in proper coordinates. Values of e < result 
in reduction in clustering strength with time. The case of 
e « — 1.2 corresponds to clusters of fixed size in co-moving 
coordinates. 

If the redshift distribution dN/dz of galaxies is known, 
then the angular correlation function w{9) can be related 
to ^(r) by the Limber equation (Limber 1953 ; Phillipps 
et al. 1978) 
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(0)=Cri"^0-('^-i)/„°°i5(z)i-^5W-i 



re 
^ Jo dz 



dz 



o-dN^^ 



Here, D(z), g{z) and C are defined as 



Diz) = 



9iz) 



c qoZ + {qo - 1)(V1 + 2goZ - 1) 



Hoq^o 



(l + z)2 



i/o(l + z)Vl + 2go^' 



C = v^. 



r ((7-i)/2) 
r(7/2) ■ 



(13) 



(14) 



(15) 



(16) 



It is evident that the exponent 7 of the spatial correlation 
function is related to the index S of the angular correlation 
function via the relation j — d + I. 

Unlike optically selected galaxies, the redshift informa- 
tion of the faint radio population is sparse and restricted 
only to sources with optical counterparts that are rela- 
tively bright. Therefore, to predict the amplitudes of the 
correlation function using equation (O) for different input 
parameters (e, Tq, 7), we need to use a model radio lumi- 
nosity function (RLF) to estimate their redshift distribu- 
tion. In the present work the evolving radio luminosity 
functions of Dunlop & Peacock (1990) at 2.7 GHz (Model 
1) and Rowan-Robinson et al. (1993) at 1.4 GHz (Model 
2), extrapolated to faint flux densities, are used. 

Dunlop & Peacock (1990) considered the free form 
luminosity function of both flat and steep spectrum ra- 
dio galaxies. Their proposed evolutionary model implies 
that the radio sources of higher power evolve faster, with 
the evolution of flat and steep spectrum radio sources 
being similar. Here we employ their RLFl model (using 
their RLF2, RLF3 and RLF4 models give similar results) 
shifted to 1.4 GHz and converted to WHz~^. The trans- 
formation from 2.7 GHz to 1.4 GHz was carried out as- 
suming a power law spectral index a {Si, ex z^~") of 0.8 
for steep spectrum and 0.0 for flat spectrum radio sources. 
The redshift distributions at different flux limits for Dun- 
lop & Peacock models are shown in Figure 0. The median 
of the redshift distribution for flux limits of 0.1, 0.5 and 
1 mJy is 0.84, 0.85 and 0.91 respectively. 

Rowan-Robinson et al. (1993) divided the radio pop- 
ulation into spirals, which they showed to be indistin- 
guishable from the starburst galaxies detected at 60 ^m 
by IRAS, and ellipticals. For the former population we 
employ Saunders et al. (1990) 'warm' IRAS galaxy com- 
ponent luminosity function at 60 /xm translated to 1.4 GHz 
using the empirical relation derived by Helou et al. (1985). 
For ellipticals, following Rowan- Robinson et al. (1993), we 




0.01 0.1 



Fig. 8. Redshift distributions for different flux density 
cutoffs using the radio luminosity function introduced 
by Dunlop & and Peacock (1990). Dashed line: Steep 
spectrum radio sources; Dotted line: Flat spectrum radio 
sources ; Solid line: both populations. The normalisation 
of the redshift distribution is arbitrary. 

use Dunlop & Peacock (1990) pure luminosity evolution 
model RLF6 (parameters from their Table C3) . To match 
the source counts and the redshift distribution of the op- 
tically brighter radio sources, pure luminosity evolution 
for the spiral galaxy component is invoked as described in 
Rowan-Robinson et al. (1993) and Hopkins et al. (1998). 
The predicted redshift distributions at different flux den- 
sity limits are plotted in Figure ^ The median of the red- 
shift distribution for flux limits of 0.1, 0.5 and ImJy is 
1.34, 1.37 and 1.44 respectively, significantly larger than 
the median values predicted by the model developed by 
Dunlop & Peacock (1990). Thus, the correlation ampli- 
tudes predicted by this RLF model are lower compared to 
the previous one, for the same parameters Tq, e. 

Using equation ( [l3| ) we predict the amplitudes of 'w{9) 
for different values of e, Tq, 7 and for different radio lu- 
minosity functions. The value of 7 is fixed to either 1.8 
or 2.1, with the models calculated for each value of 7 in- 
dependently. Upper limits for the clustering length Vo for 
e = —1.2, are listed in Table 0. The error estimates of Tq 
are calculated from the •iapoisson correlation amplitude 
uncertainties (Table H). 

7. Discussion 

Peacock & Nicholson (1991) have shown that apparently 
bright radio sources {S > 0.5 Jy) in the redshift range 
0.01 < z < 0.1 are spatially correlated, following a power- 
law behaviour ^ = (r/llh~^ Mpc)~^-^. Loan, Wall & La- 
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Fig. 9. Redshift distributions for different flux density 
cutoffs using the radio luminosity function described by 
Rowan-Robinson et al. (1993). Dashed line: Ellipticals; 
Dotted line: Spirals; Solid line: both populations. The nor- 
malization of the redshift distribution is arbitrary. 
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Fig. 10. Predicted amplitude of the angular correlation 
function at different flux density cutoffs, assuming that 
the radio population consists of starbursts and ellipti- 
cal galaxies with different clustering properties Solid line: 
To = 5h~'^Mpc for starbursts, Tq = llh~^Mpc for el- 
lipticals (Model A); Dashed line: ro=llh^^ Mpc for both 
populations (Model B). The clustering evolution param- 
eter is taken to be e=-1.2 and 7=1.8 for both popula- 
tions. Also shown are the estimated amplitudes from this 
work (filled circles; upper limits) and for sources with 
1 < S'i.4 < 2mJy detected in the FIRST survey. 



Table 4. Upper limits of the clustering length, for e=- 
1.2 and for two RLF models. Model 1: Dunlop & Peacock 
(1990) RLF; model 2: RLF described by Rowan-Robinson 
et al. (1993). The superscripts are the upper limits to the 
errors in the Tq estimates, corresponding to the bootstrap 
uncertainties of the angular correlation amplitude. 



hav (1997), using two large-area 4.85 GHz radio surveys 
covering 70% of the sky, conclude that the 2-D distribution 
of radio sources brighter than 50 mJy is consistent with a 
correlation length in the range 13 < r,, < 18 h^^Mpc 
and an evolution parameter —1.2 < e < 0. This result is 
based on the RLF models developed by Dunlop & Peacock 
(1990) and a value of 7=1.8. In a separate study. Cress 
et al. (1996) estimated the angular correlation function 
of the FIRST radio survey for sources with flux density 
5*1.4 >lmjy. Assuming e = —1.0 and 7 = 2.1, they found 
r„=6-8h~iMpc (Cress et al. 1997). 

The To upper limits estimated in this study at different 
flux density cutoffs agree with the estimates derived from 
brighter radio samples (Cress et al. 1997; Loan, Wall & 
Lahav 1997). Because of the small number of sources in 



each of the independent flux-limited sub-samples («318) 
the uncertainties are large and thus we cannot conclude 
from the present sample if there is a change in the clus- 
tering properties of radio sources with flux density. Figure 
EO compares the correlation amplitude upper limits from 
this study with those calculated for the radio sources de- 
tected in the FIRST radio survey with 1 < 5'i.4 < 2mJy 
(Cress et al. 1996). The agreement is good, even within 
the Poisson errors. 

In recent years, deep radio surveys, (e.g. Windhorst et 
al. 1985; Fomalont et al. 1997; Hopkins et al. 1998) have 
shown a flattening in the slope of the normalised number 
counts at sub-mJy levels, revealing an excess of faint radio 
sources over the 'normal' radio population of giant ellipti- 
cals and QSOs. To address this problem, models invoking 
strong evolution of either spiral galaxies (Condon 1989) or 
star-forming IRAS galaxies (Rowan-Robinson et al. 1993; 
Hopkins et al. 1998) have been combined with RLF mod- 
els of the 'normal' radio population. This scenario is sup- 
ported by photometric and spectroscopic studies reveal- 
ing that the sub-mJy radio sources can be identifled with 
galaxies exhibiting evidence of increased star formation 
(Benn et al. 1993; Windhorst et al. 1985). Furthermore, 
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studies of local galaxies with enhanced star formation (late 
spirals, IRAS galaxies, HII galaxies), has shown that those 
objects are more weakly clustered {ro ss2-4h^^Mpc, for 
7 « 1.5 - 1.7; Davis & Geller 1976; Giovanelh et al. 1986; 
Saunders, Rowan-Robinson & Lawrence 1992 ) than E/SO 
galaxies. This implies that the sub-mJy population could 
be more weakly clustered than the E/SO objects that host 
the majority of apparently brighter radio sources. Support 
for this view was advanced by Cress et al. (1996) who 
found that the slope of the angular correlation function 
of the sources with 1< S'1.4 < 2mJy is flatter compared 
to that found for S'1.4 >3mJy. They interpret this result 
as being due to the increased contribution from starburst 
galaxies at lower flux density limits, which have flatter an- 
gular correlation functions compared to ellipticals, which 
dominate at brighter flux densities. A similar argument is 
used by Peacock (1997) to explain the apparent conflict 
between the value of Tq =6.5h^^Mpc found for a sam- 
ple of radio sources brighter than 2.5 mJy and the value 
To =11 h~^ Mpc found by Peacock & Nicholson (1991) for 
sources brighter than 500 mJy and redshifts in the range 
0.01 < z <0.1. 

To explore further the implications of the two scales 
of clustering, and in particular to explore the potential 
to eliminate competing models by observations, the am- 
plitude of the angular correlation function was estimated, 
by adopting a simple model (Model A) in which the radio 
population consists of two components, one dominating 
at brighter (>lmjy) fluxes, with correlation length Vq — 
11 h^^ Mpc (7=1.8) and the other dominating at sub-mJy 
levels with Tq = 5h~^ Mpc (7=1.8) similar to that found 
for local starburst galaxies. Any cross-correlation between 
the two radio populations is ignored. The clustering evolu- 
tion parameter is taken to be e=-1.2 and our RLE model 
2 is employed to predict the redshift distribution of the 
two radio populations at faint flux densities. In Figure |o| 
we plot the flux density cutoff against the amplitude of 
the angular correlation function calculated from model A. 
Eor comparison, the expected relation assuming the same 
value of To = 11 h~^ Mpc (7=1.8) for the two populations 
(Model B), is also plotted, along with the upper limits for 
the angular correlation amplitudes, A^,, calculated in the 
present study for different flux density cutoffs. The uncer- 
tainties are too large to allow discrimination between the 
two models. 

It is an interesting exercise to predict the depth and 
the solid angle subtended by a radio survey that would re- 
veal at a 3(7 significance level if the clustering properties of 
faint radio sources were significantly different from those 
of the brighter ones as a result of the changing popula- 
tion. This is accomplished by estimating the uncertainty 
in A^ for a survey of a given solid angle and completeness 
limit, as described in Appendix A. This then is compared 
to the difference Aj4u, between the correlation amplitudes 
predicted at the same flux density cutoff from Models A 
and B. To discriminate between the two models at a 3cr 





^ 


1 1 


^ 




1 


MM 


1 1 1 1 1 1 M 1 II 




; 












^0.5"^*^ ; 


^ 1 














^ '' 





- 












^z'' ^deg — 


6 


= 














- 












^^^ - 




' 












/'' y'' 




" 












/ ^-^ 




- 












^3dcg _ 



































d 
















: 








7^ 


~-~~^^^ 


^ ^' : 




















I 












^-\.'' IQdeg - 




- 






,■-' 




■■" 


,''\^ -'''' 


^ 


- 


^ 


, '' 























'' 


•;><C 





' 

















I 










■' 


_,■' ^^ ~^-^~^ - 


6 


- 












_,'■ ^-^ 




















■ ' 












X -^ 




' 












/' 


tn 




^i 


iIa 




1 ^ 


1 1 1 n1 


\ 1 1 1 1 1 1 1 1 II 



0^01 0^1 1 

S (mJy) 

Fig. 11. Flux density cutoff versus the estimated uncer- 
tainty in the amplitude of the angular correlation function. 
Dashed lines correspond to surveys of circular geometry 
with different radii. The Phoenix survey has a radius of 
1°, while the FIRST radio survey is subtending, at the 
moment, an area corresponding to a radius of r::22°. The 
solid line represents the locus of 3cr confldence level dis- 
crimination between models A and B described in the 
text. 



confldence level, the uncertainty in A^ should be 3 times 
smaller than AA^,. The results are shown in Figure ^ 
where the flux density cutoff is plotted against the un- 
certainty in Aw for different survey areas. The solid line 
delimits the area in the parameter space in which bA^ 
is at least 3 times smaller than A^^; and hence deflnes 
the locus of 3(T confidence level discrimination between 
the simplified models A and B. We conclude that either a 
very deep survey (5;im=40/iJy) over « 1 deg^ or a survey 
over a larger area at a brighter limit (e.g. at S'/im=0.2 mJy 
over « 27deg^) is required to discriminate between models 
A and B and thus reveal if there exists a weakly clustered 
radio population at faint fiux densities. Our analysis shows 
that the FIRST radio survey, (Sum — 1 niJy) covering, at 
the moment, an area of « 1500 deg^, is also appropriate 
for this purpose. 



8. Conclusions 

Using a deep, homogeneous survey covering an area of 
~ 3deg^, we have investigated the clustering properties of 
faint (sub-mJy) radio sources by calculating the angular 
correlation function of the sample. Simulations have been 
carried out to investigate the significance of our results. 
It has been demonstrated that the amplitudes estimated 
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here imply a non-uniform distribution of radio sources, 
albeit at the 2cr significance level. 

Furthermore, we find that the correlation amplitude 
upper limits estimated here, despite the large uncertain- 
ties, are consistent with the results derived in other stud- 
ies from larger samples at brighter flux density limits. 
Spectroscopic and photometric studies of sub-mJy radio 
sources show that the majority have properties similar to 
starburst galaxies. Since local galaxies, with increased star 
formation, exhibit weaker clustering than average, it is 
tempting to link the clustering properties of the sub-mJy 
population with those of local starburst galaxies. However, 
the present data, due to small number statistics, cannot 
reveal a change in radio source clustering properties with 
flux density. Nevertheless, we have demonstrated that fu- 
ture surveys having sensitivities at fiJy flux density levels 
and covering large areas on the sky, will be suitable to ex- 
plore possible changes of the correlation amplitude with 
flux density. 
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Appendix A: Angular Correlation Function 
Amplitude Uncertainty 

The uncertainty in the amplitude of the angular correla- 
tion function of a fiducial radio survey with limiting flux 
density Sum, subtended over a solid angle fi, is estimated 
here. To simplify the procedure, the following assumptions 
are made: 
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1. the area subtended by the survey has a circular geom- 
etry; 

2. the angular correlation function of radio sources is a 
power law of the form w{9) = Au,0^^ , with the index 
6 being fixed; 

3. the values of w{9) at different angular separations are 
uncorrelated; 

4. The uncertainty in w{6) follows Poisson statistics (Pee- 
bles 1980) and is given by 



M^)--'' + "^'^ 



DD ' 

(Al) 



where DD is the number of data-data pairs with an- 
gular separations 6 and 6 + d6, N is the number of 
sources above the flux density limit and < dQ > is the 
mean value of solid angle subtended by annuli with 
inner and outer radii 9 and 9 + d9 respectively. The 
value oi < dfi > is affected by the field boundaries. In 
the present study, this is estimated by creating simu- 
lated catalogues of sources distributed randomly over 
the area of the survey and calculating the expected 
number of random points, RR, with separations 9 to 
9 + d9. This is given by 

RR=^NniNn-l)^^, (A2) 

noindent where Nn is the number of sources in 
the simulated dataset.The assumption of Poisson er- 
rors underestimates the uncertainty in A^, (Bernstein 
1994). 

The assumptions (iii) and (iv) allow us to use least 
squares fitting techniques to estimate the uncertainty in 

-^(,£(fe)T^ 

where wn is the integral constraint defined in section y. 
The summation is carried out between a minimum {9min) 
and a maximum {9 max) angular separation within equally 
separated logarithmic bins of width Alog6' = 0.2 [9 in 
degrees). The former is taken to be 3' and the latter is set 
three times smaller than the radius of the survey, while it 
cannot be larger than 1°. This is because in real surveys 
one needs to minimise the effects of field boundaries when 
estimating A^ . Additionally, one is interested in the value 
of All, on small scales, while for larger scales the angular 
correlation function may deviate from a power law. For 
example in the optical wavelengths Maddox et al. (1990a) 
found a break in w{9) at scales 9 « 1.5". 



